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Abstract The Phoenix Deep Survey (PDS) is a multiwavelength survey based on deep 
1 .4 GHz radio observations used to identify a large sample of star forming galax- 
ies to z = 1. Here we present an exploration of the evolutionary constraints on 
the star-forming population imposed by the 1.4 GHz source counts, followed by 
an analysis of the average properties of extremely red galaxies in the PDS, by 
using the "stacking" technique. 
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1. Introduction 

The study of galaxy evolution in recent years has included a strong focus on 
the star formation properties of galaxies. Many of these studies are based pri- 
marily on selection at ultraviolet (UV) and optical wavelengths, and are known 
to be strongly affected by obscuration due to dust. It has been shown that selec- 
tion at these wavelengths results in samples of star forming systems that omit 
a significant fraction of heavily obscured galaxies [11]. There have in addition 
been suggestions that the most strongly star forming systems suffer the most 
obscuration [1, 8, 3, 12, 7]. Using radio selection to construct a star forming 
galaxy sample allows the detection of such systems, and the average obscura- 
tion in this case is significantly higher than in optically selected samples [1,8]. 

To identify a homogeneously selected catalogue of star forming galaxies, 
unbiased by obscuration effects, and spanning a broad redshift range (0 < z < 
1), the Phoenix Deep Survey (PDS) uses a deep 1.4 GHz mosaic image, based 
on observations with the Australia Telescope Compact Array. This has been 
used to construct one of the largest existing deep 1 .4 GHz source catalogues [6] 
(see http://www.atnf.csiro.au/people/ahopkins/phoenix/) from which the star 
forming galaxy sample will be drawn. The PDS has already been highly sue- 
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Figure 1. Likelihood contours in the (P, Q) plane, indicating the degeneracy between rates 
of luminosity and density evolution for star-forming galaxies. 



cessful in providing a basis for several investigations of the nature of star form- 
ing galaxies and their evolution (see references in [6]). Throughout the present 
investigation we assume a (Qm = 0.3, = 0.7, Hq = 70) cosmology. 

2. 1.4 GHz source counts and models 

The PDS 1.4 GHz source counts have already been explored in detail [6]. 
Here we explore source count models based on local luminosity functions 
(LFs) and their assumed forms of evolution. We assume an evolving LF for 
AGNs [5] (converted to our assumed cosmology), since we are interested here 
in exploring the properties of the star-forming (SF) galaxies. We then use a 
measured local 1.4 GHz LF for SF galaxies [10], and compare the observed 
source counts with the model predictions subject to a range of both luminosity 
evolution, L(z) oc (1 + z)®, and density evolution, <f>(z) oc (1 + z) p , for the 
SF population. The \ 2 estimator for each combination of (P, Q), gives the sta- 
tistical likelihood, and contours showing the region of maximum likelihood in 
the (P, Q) plane are shown in Figure 1 . There is a clear degeneracy, with max- 
imum likelihood occuring for any {P,Q) satisfying Q s=s 2.7 — 0.6 P (note that 
only positive (P, Q) have been considered in this analysis). This sort of degen- 
eracy has also been shown in investigations of optically selected galaxies [9]. 
To illustrate the effect of the degeneracy, Figure 2 shows the predicted source 
counts compared with the observations for five (P, Q) pairs, indicated by the 
positions a-e in Figure 1. The redshift distributions predicted by these models, 
however, are significantly different, and photometric redshifts are now being 
utilised in order to differentiate between the different evolutionary models. 
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Figure 2. Observed 1.4 GHz source counts compared with model source count predictions 
(solid lines). The models, while all predicting similar source count distributions, span a broad 
range of combinations for the luminosity and density evolution parameters. 



3. Extremely red galaxies 

A sample of over 400 extremely red galaxies (ERGs) has been compiled 
from PDS observations, using the colour criterion R — K > 5. Of these, 
about 20 are detected at 1.4 GHz, while the majority (over 90%) remain un- 
detected. In an effort to understand more about the properties of these sys- 
tems, the "stacking analysis" technique often used at X-ray wavelengths has 
been adopted. We extract sub-images from the radio mosaic at the location of 
the non-detected ERGs, and construct the weighted average of the sub-images 
(weighted by 1/rms 2 , to maximise the resulting signal-to-noise, since the radio 
mosaic has a varying noise level over the image). Sub-images where low S/N 
emission (> 1.5 a) is present at the location of the non-detected source are 
excluded from the stacking, in order to avoid biasing the stacking signal result 
by the presence of a small number of low S/N sources. A further 22 (out of 
399 candidate non-detected ERGs) were excluded in this manner, leaving 377 
to contribute to the stacking signal. The stacked image, shown in Figure 3, 
has an rms noise of 1.1 /uJy, and a ps 6 a detection at 6.5 fiJy. While the ac- 
tual redshift distribution of these objects is unknown, a significant fraction of 
ERGs appear to be dusty starbursts at z > 1 [11]. Assuming that the average 
redshift of these sources is z ps 1, the inferred average ERG 1.4 GHz luminos- 
ity is 3.5 x 10 22 WHz -1 . This corresponds to an average star formation rate 
(assuming the ERGs are all star-forming systems) of ps 2OM yr~ 1 (adopting 
the recent calibration of Bell 2003 [2]; this value doubles if the calibration 
of Condon 1992 [4] is used). Clearly these results rely on a large number of 
assumptions, but they serve the useful purpose of providing a preliminary esti- 



Figure 3. Stacked image, incorporating 377 ERO candidates, from the PDS radio mosaic. 
The image rms is about 1.1 /ijy, and the stacked signal is detected at « 6 a with a flux density 
of 6.5 /j,Jy. 

mate for the average properties of these systems in order to support proposals 
for further, more detailed investigation. 
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